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Radio  propagation 
Waveguide  modes 
Low-frequency  propagation 


<^A8STRACT  /Continue  on  leve'se  -t  necessmy  Opd  identity  by  block  number) 

—‘For  low-frequency  (W)  cases  of  nighttime  propagation,  an  elevated  antenna,  or  distances  close  to  the  transmitter, 
the  propagation  constants  for  a  large  number  of  modes  must  be  found  to  adequately  represent  the  fields  in  the  earth- 
ionosphere  waveguide.  Most  of  the  cost  in  using  MODESRCH  to  obtain  these  mode  constants  has  been  in  performing 
full-wave  integrations  to  obtain  the  reflection  matrix,  Until  recently  it  had  been  necessary  to  carry  out  several  full- 
wave  integrations  each  for  and  its  derivative  with  respect  to  the  angle  of  incidence,  6,  for  each  waveguide  mode.  As  a 
result  of  this  present  effort  it  has  become  possible  to  obtain  interpolated  values  of  jk  which  are  more  than  adequate  to 
describe  the  earth-ionosphere  mode  constants. 

The  approach  used^  is  to  transform  the  reflection  matrix,  R,  obtained  from  a  full-wave  integration,  into  magnetoionic 
reflection  coefficients,  .  Each  coefficient  represents  the  complex  amplitude  of  the  ordinary  or  extraordinary  component 
of  the  reflected  wave  given  an  incident  wave  which  is  polarised  at  either  ordinary  or  extraordinary.  The  equations 
defining  the  magnetoionic  eigennvectors,  which  are  required  for  this  transformation,  have  been  derived  for  the  free-apace 
limit,  so  that  they  are  applicable  to  a  height  that  it  below  most  of  the  ionisation. 

Experience  infuse  of  the  algorithm  has  shown  that  for  eigenangles  not  too  near  branch  points,  the  complex  logs  of 
the  elements  of  R  vary  so  nearly  linearly  that  values  of  R  derived  from  interpolated  values  of  fnR  *  are  more  than 
adequate  for  use  In  determining  mode  constants.  An  example  is  given  for  illustration,  (continued  next  page) 
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low>frequency  cases  requiring  a  large  number  of  modes  results  in  considerable  cost  savings. 
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I.  INTRODUCTION 


For  low-frequency  (LF)  cases  of  nighttime  propagation,  an  elevated  antenna,  or 
distances  close  to  the  transmitter,  the  propagation  constants  for  a  large  number  of  modes 
must  be  found  to  adequately  represent  the  Fields  in  the  earth-ionosphere  waveguide.  Most 
of  the  cost  in  using  MODESRCH  (Morfitt  and  Shellman,  1976)  to  obtain  these  mode  con¬ 
stants  has  been  in  performing  full-wave  integrations  to  obtain  the  reflection  matrix,  R.  Until 
recently  it  was  necessary  to  carry  out  several  full-wave  integrations  for  R  and  its  deri¬ 
vative  with  respect  to  the  angle  of  incidence,  8,  for  each  waveguide  mode. 

A  means  of  using  interpolated  values  of  R  is  obtained  in  MODESRCH.  The  full- 
wave  solutions  are  found  at  the  four  corners  of  each  of  a  number  of  rectangles,  and  approx¬ 
imate  values  of  elements  of  R  are  found  using  a  third-order  interpolation.  The  functional 
forms  of  the  elements  of  R  are  sufficiently  complicated,  however,  that  they  cannot  be 
adequately  approximated  by  the  interpolation  except  for  use  in  finding  preliminary  values 
used  internally  in  MODESRCH. 

It  was  perceived  that  some  of  the  nonlinear  variation  of  elements  of  R,  especially 
in  the  LF  range,  might  be  due  to  interference  of  waves  reflected  from  two  fairly  distinct 
complex  altitudes.  The  relative  phases  of  reflections  from  the  two  altitudes  depend  strongly 
on  the  angle  of  incidence,  6.  The  approach  taken  in  the  formulation  described  in  sections  II 
and  III  is  to  separate  the  reflection  coefficients  into  ordinary  and  extraordinary  components 
at  a  height  that  is  below  most  of  the  ionization.  Since  neither  component  is  expected  to 
vanish,  the  complex  logs  of  the  elements  of  the  magnetoionic  reflection  matrix  may  be  used 
in  the  interpolation. 

At  each  branch  point  there  is  an  ambiguity  in  separating  the  ordinary  and  extra¬ 
ordinary  components.  A  formulation  for  locating  these  points  in  the  0-plane  is  described 
in  section  IV.  Near  a  branch  point  the  elements  of  R  are  used  as  in  MODESRCH. 

There  is  an  ambiguity  in  the  relative  number  of  whole  cycles  in  the  complex  logs  of 
the  magnetoionic  reflection  coefficients  at  the  four  corners  of  the  search  rectangle,  (figure  1 ). 


Figure  1.  Search  area. 
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This  is  resolved  by  choosing  the  numbers  of  cycles  so  as  to  minimize  the  magnitudes  of  the 
second*  and  third-order  coefficients  of  the  interpolation  series.  The  formulation  for  obtain¬ 
ing  the  coefficients  and  for  making  the  test  is  described  in  section  V. 

An  example  of  use  of  the  formulation  presented  in  this  report  is  given  in  section  VII. 
The  frequency  used  is  60  kHz,  and  a  nighttime  ionosphere  model  is  used. 


II.  MAGNETOIONIC  COMPONENTS 

Decomposition  of  the  reflection  coefficients  into  ordinary  and  extraordinary  com¬ 
ponents  requires  the  magnetoionic  eigenvectors  at  a  height,  z^,  that  is  below  most  of  the 
ionization,  that  is,  in  the  limit  of  small  values  of  electron  density.  Since  an  earth  curvature 
term  is  included  in  the  full-wave  solution  for  R,  this  term  is  also  included  in  the  form  for 
the  eigenvectors. 

The  differential  equation  for  the  propagation  of  radio  waves  in  the  ionosphere,  in 
the  vertical  direction,  is  given  by  Budden  (1961)  as 

?'  =  -i«T?  ,  (1) 

where  the  elements  of?'  are  Ex.  -F.y.  5fx.  3fy.  The  matrix,  T.  is  defined  by  Budden  (1961) 
and 

i  =  (-1  Y1 
k  =  wave  number. 

The  prime  denotes  differentiation  with  respect  to  height,  z. 

Equation  ( 1  »  may  also  be  written  as 

?’  =  -ikT  ?  . 

where  the  elements  of  w  are  labeled,  in  order.  7fy.  Ey,  3fy  and  Ey.  Also, 

T  =  L_l  T  L  . 

vv  here 


and 

q  =  (C-  +  el’ :  Re  (q  I  '  *  0 
(  -  co ,  0 


X 


h  =  reference  height  for  earth  curvature 
re  =  radius  of  the  earth 

=  height  at  which  magnetoionic  components  are  evaluated 

The  matrix,  L,  is  defined  in  a  way  analogous  to  that  of  Budden  (1961)  but  such  that  the 
superscripts  u  and  d  refer  to  upgoing  and  downgoing  waves,  respectively,  at  the  height 
z  =  zj,  where,  in  general,  e  =£  0.  The  branch  cut  for  q  is  shown  in  figure  2. 


Figure  2a.  Branch  cut  for  zd  >  h. 


Figure  2b.  Branch  cut  for  zd  <  h. 


At  this  height  and  with  a  vanishingly  small  value  of  electron  density,  the  matrix  T 
approaches  being  a  diagonal  matrix  with  elements  q,  q,  -q,  and  -q.  Hence,  the  characteris¬ 
tic  upgoing  waves  at  z  =  z^  are  uncoupled  from  the  characteristic  downgoing  waves.  How¬ 
ever.  since  for  either  upgoing  or  downgoing  waves  the  two  eigenvalues  of  T  are  the  same 
(equal  to  q  for  upgoing  waves  and  to  -q  for  downgoing  waves),  it  is  necessary  to  consider 
the  terms  of  T  which  are  present  for  a  very  small,  but  otherwise  negligible,  value  of  electron 
density  in  order  to  determine  the  polarization  of  characteristic  waves  at  the  bottom  of  the 
ionosphere  model.  For  this  purpose  only  the  upper  left  2X2  submatrix  of  T  and  the  lower 
right  2X2  submatrix  of  T  are  needed.  These  are  used  in  the  pair  of  matrix  equations 


w!.h  =  -i*c  T ,„.w,lH 


where  the  subscript  ud  refers  to  there  being  two  sets  of  variables,  one  for  upgoing  waves 
and  one  for  downgoing  waves. 

To  form  the  matrices  TU(j,  the  susceptibility  matrix  for  very  small  values  of  electron 
density  is  needed.  This  may  be  written 


M  =  C'« 


-y-G-i/G 


-n-tmG 


m-fcnG 


n-CmG 


-m  ~G- 1  /( i 


-C-mnG 


-in-CnG 


fe-mnG 


-n-G-l/G 


where 


(  M  -  iX/[  U“Y(  U~  -  Y~)l 
G  =  Y/(Z  +  it 

X  =  normalized  electron  density 
Y  =  normalized  magnetic  field  strength 
Z  =  v‘u) 

v  =  collision  frequence 
co  =  angular  wave  frec|uency 
l'  -  1  -  iZ 


fe  =  sin  6  cos  cv 


m  =  sin  S  sin  a 


n  =  -  cos  6 


h  ■  codip  angl  • 

o  -  azimuth  ol  propagation  measured  east  of  north 


Note  that  the  value  of  v  must  correspond  to  the  height,  z^.  which  is  below  most  of  the  ioni¬ 
zation.  1  or  MODI  SRCII.  the  reference  height.  It.  for  earth  curvature  has  usually  been  taken 
to  lie  SO  km.  but  it  ma\  be  set  to  the  height  of  the  base  of  the  ionosphere  or  to  some  other 
heinht .  I  hen 


Si  m  -t  nf<  i )  p~ 


S(  f  -  mnG 


-( n  -  VmG ) 


Q  -  nr(i  -  |/G 


-StC  +  mnGl/p 


r  -  v-t, -  i  (> 


it  +  tiiiti 


-S(m  -  nt!G)  p 


14  -  <U  -  n~G  -  1  G )  p-  +  c|“t  n“ G  +  I  /( i ) /p“ 


Finally, 


r  -  p2/Cm 

Q  -  i2/cm- 


t-bzg 


T  =  Cq 


A  -  mBG 


where 


CT  =  ±udCM/(:<^P‘') 
A  =  CS±udnq 
B  =  -nS±utj8q 

T  =p2(2Q-l/G) 


±ud  =  +  ^or  uPS°'ng  waves 
=  -  for  downgoing  waves 


The  condition  for  determining  eigenvalues  and  vectors  is  then 
i2r._-v  _tA  +  mRran2'  w  '  o' 


-B~G-X  -(A  +  mBG)pz  !Cy 


q  q 

A  -  mBG  -m"p-G  -  X 


The  eigenvalues  are  found  from  a  solution  of  the  characteristic  equation 


X- +  (B_  +  m-p“)GX  +  A_p“  =  0 


Using  the  first  row  of  the  matrix  in  (2).  the  eigenvectors  at  the  base  ot  the  iono¬ 
sphere  for  either  upgoing  or  downgoing  waves  are  then  given  by 


:( 1 )  — 


:ud 


where 


3fy  j  -  2(A  +  mBG)p2 
1  yl  =w  Vf 


and  where 


W  =  (m~p“  -  B“)G 

f  =[(B2  +  m2p2)2G2-4A2p2]‘/2  •  (3) 

It  is  the  choice  of  the  sign  ±oe  associated  with  f  that  distinguishes  the  “ordinary”  waves 
from  the  "extraordinary”  waves. 

Using  the  second  row  of  the  matrix  in  (2),  the  matrix  of  eigenvectors  may  also 
be  written 


where 

*y2  =  w  +oer 

1  v->  =  2(  A  -  mBG> 

Neither  E^j*  nor  is  in  normalized  form.  However,  it  can  be  seen  that  the  two 
column  vectors  are  given  in  tire  same  order  in  the  two  matrices  by  noting  that 

r2  =  W:-3Cyl  Ey  ,. 

The  form  that  is  most  manageable  is  derived  by  using  the  left  column  vector  of  and  the 
right  column  vector  of  f  j and  normalizing  the  matrix  so  that  the  determinant  is  equal  to 
unity.  This  results  in 


where  the  sign  of  the  square  root,  representing  the  value  of  f  in  (3).  is  chosen  so  that  the 
sum  W  +  f  has  the  larger  magnitude  in  any  one  search  rectangle.  The  eigenvector  matrix, 
Euj.  is  singular  for  values  oft?  for  which  f  =  0,  that  is.  for  which  ordinary  waves  cannot 
be  distinguished  from  extraordinary  waves.  Although  not  singular,  it  is  not  well  conditioned 
for  large  collision  frequency  when  the  propagation  is  east-west  at  the  equator. 

Note  that  for  a  small  value  of  G.  that  is.  for  a  large  value  of  collision  frequency  at 
the  height  z(j  but  not  for  east-west  propagation  at  the  equator,  the  matrix  of  eigenvectors 
is  approximated  by 


6 


£ud  (large  v)  ^ 


s/2  . 


representing  nearly  circularly  polarized  characteristic  waves. 

For  east-west  propagation  at  the  equator.  £  =  n  =  0.  In  this  case  EU[j  becomes 


Eud  (equator  east-west) 


4 

l  hc  let t  eigenvector  corresponds  to  ordinary  waves  since,  for  this  case.  5Cy  =  0.  For  extra¬ 
ordinary  waves  l -  =  0  for  east-west  propagation  at  the  equator  ( Budden,  1961). 


Ill  REFLECTION  COEFFICIENTS 

Once  the  matrices  of  eigenvectors  Eu  and  Ej  are  known,  the  elements  of  Roe  may 
be  found  from  those  of  R  and  vice  versa.  The  reflection  matrix.  R.  is  as  defined  by  Budden 
i  19(>1  i.  It  the  first  columns  of  E  and  Ej  correspond  to  ordinary  waves,  then  Roe  is  of  the 
form 


where  the  first  subscript  refers  to  the  polarization  of  the  upgoing  waves  and  the  second 
subscript  refers  to  that  of  the  downgoing  waves. 

Matrices  C<K>  and  I)oc  are  used  to  define  upgoing  and  downgoing  waves,  respec¬ 
tively.  where  the  elements  of  the  first  row  of  each  represent  values  of  and  the  elements 
of  the  second  row  of  each  represent  values  of  I  y  The  first  columns  of  U°°  and  Doe  cor¬ 
respond  to  the  condition  that,  for  the  upgoing  wave.  =  1  and  Ev  =  ().  and  the  second 
columns  correspond  to  the  condition  that,  tor  the  upgoing  waves.  7fy  =  0  and  Ev  =  1.  Then 


Doe  =  E".1  R 

-■'w  -w 

and 

Roe  =  |)Oe  (Uoerl  _  E-1  rj 

^  ^  'vvll 


The  values  are  used  in  the  interpolation  scheme  as  “given"  values  at  the  upper  corners  of 
each  search  rectangle. 


Approximate  values  of  the  elements  of  Roe  found  by  interpolation  must  be  trans¬ 
formed  to  values  of  elements  of  R.  The  first  step  is  to  form 

y  =  £u 

5  =  Ed  Roe. 

The  elements  of  the  first  rows  of  U  and  D  represent  values  of  Hy,  and  the  elements  of  the 
second  row  of  each  represent  values  of  Ey.  The  first  columns  of  U  and  D  correspond  to  the 
condition  that  the  upgoing  wave  is  ordinary,  and  the  second  columns  of  U  and  D  correspond 
to  the  condition  that  the  upgoing  wave  is  extraordinary,  if  the  first  columns  of  Eu  and  E^ 
correspond  to  ordinary  waves.  Then 

R  =  DU-1  =  Ed  Roe  H~' . 

Two  effective  complex  heights  of  reflection  may  be  found  by  considering  that  the 
diagonal  elements,  rk.  of  Roc  would  be  expected  to  vary  as 

drk/dC=  2i«(lik-/d>rk  . 

where  rj  -  ROQ.  r->  =  R  ,  zd  is  the  height  at  which  Roe  is  defined,  and  x  is  the  wave  num¬ 
ber.  Then  the  effective  heights  of  reflection  are  given  by 

hk  =  zd  +  {d[Cn(rk)!/dC}/(2ix) 

It  would  be  expected  that  the  ordinary  wave  is  reflected  from  the  upper  height  (Budden, 
ldbl  )  and  that  the  reflection  from  this  height  is  the  stronger  of  the  two. 


IV.  BRANCH  POINTS 

Hie  reflection  matrix  R  cannot  be  resolved  into  ordinary  and  extraordinary  com¬ 
ponents  at  points  in  the  complex  0-plane  where  the  argument  of  the  square  root,  which  is 
the  expression  for  f.  vanishes.  Near  these  points  the  reflection  matrix  R  must  be  used  in  the 
search  for  the  waveguide  eigenangles  without  a  transformation  to  the  Roe  form. 

I  he  expression  (3)  for  f  in  section  11  may  be  written 


r  =  o+i->': 


w  here 


f±  -  (B-+m*-p“Ki  t  2Ap  . 

and  where  Ret  p)  >  0.  1  here  is  a  branch  point  at  each  value  of  6  for  which  either  f+  or  f 
vanishes.  This  pair  of  conditions  vields 

,  ^  3 


1 .  =  |  p--(  f  Sj  +  nq j  1“  |  ( ■  *  2p(  fSj  +  in],  I  -  U 


(4) 


m  m  nr w  mx  HA  M  W  iff  iff  iff  iff  iff  iff  iff  iff  iff  iff  iff iff  iff  iff  w*  in  g*  iff  iff  jui  i*  iff  i»  iff  iff  w' ■ 


for  which  a  partial  solution  is 

(CSj  +  nqj)  =  tj  =  p  1 1.0  ±  (1  +G2)'/2]/G 

=  -p  [i.o  ±o  +g2),/2]/g 

where 

I/. 

Re(Sj)  >  0  . 


i  =  1,2 
i  =  3,4  , 


Sj  =  (p--qf)/2 


Note  that  the  ±U(j  sign  has  been  omitted  so  that  values  of  q  for  both  upgoing  and  downgo¬ 
ing  waves  are  included  in  (4).  The  values  of  q  at  the  branch  points  are  then  the  solutions  to 


fj  =  tj  -  (CSj  -  nqj)  =  0  . 


The  solutions  for  C  =  0  are  q:  =  tj/n.  For  daytime  ionosphere  models.  G  is  small  at 
the  base  of  the  ionosphere  model  since  Z  »  Y.  In  the  limit  of  large  collision  frequency. 

G  -*■  0.  and  the  two  solutions  are  q:  =  0.  For  this  case.  then,  there  are  branch  points  at 
0  =  cos“'<  ±ie' :).  These  are  usually  within  about  10  deg  of  the  90-deg  point  in  the  0-plane. 

The  solutions  for  a  nonzero  given  value  of  6.  which  is  a  function  of  the  given  values 
of  codip  and  azimuth  of  propagation,  are  found  by  solving  (5)  for  successively  larger  values 
of  6.  Use  is  made  of  the  derivatives 


dfj/dqj  =  n  -  Cqj/Sj 


d“fj/dq“  =  -  (6/ Sj)  ( 1  -  qf/Sp. 


i  i 


For  each  new  value  of  C.  the  increment  in  q,  is 
Aqj  =  -fj  '(dfj/dqj) 


dfj'dqj  »  fjid-fjdqp 


-y,  -  <yf  -  -fj*  ; 


1  M 

otherw  ise 


where 


y'j  =  (dfj/dqj  )/(d“fj/dqp 


and  where  the  sign  is  chosen  that  yields  the  smallest  magnitude  of  Aqj. 

I  he  search  is  terminated  when  t  becomes  equal  to  its  originally  specified  value 
Note  that  .with  the  condition  Re(  S)  >  0,  this  yields  four  values  of  q  j  The  branch  points 

are  then  given  by 


-1 


lj  =  cos  1  Cj 


q 


where 


The  sign  of  Cj  is  chosen  so  that  the  real  part  of  Cj  has  the  same  sign  as  the  real  part  of  qj. 

If  the  real  part  is  positive,  the  branch  point  pertains  to  upgoing  waves  and  to  Eu.  Otherwise 
it  pertains  to  downgoing  waves  and  to  E^. 

V.  INTERPOLATION 

The  interpolation,  within  each  rectangle,  is  based  on  the  reflection  coefficients 
found  from  the  full-wave  solution  at  the  corners  of  the  rectangle.  The  four  sets  of  reflection 
coefficients  provide  for  interpolation  in  the  third-order  form 

R  =  R +  R'.  <0  -  0,,)  +  R"  (0  -  0.)2/2  +  R'"  (0  -  0L.)3/ 6  . 

where  0C  is  the  value  of  0  at  the  center  of  the  rectangle.  The  values  of  Rc  and  its  derivatives 
are  considered  to  be  in  terms  of  coefficients  such  that 

Sc  =  ^0 
&  = 

K  =  2as/t2 
R"'  =  6a,'t3  . 

where  t  is  a  real  number  representing  the  length  in  the  6  plane  from  the  center  of  the 
rectangle  to  any  one  of  the  corners. 

f  or  simplicity,  the  matrix  notation  is  dropped  at  this  point  and  the  scalar  variables 
are  taken  to  refer  to  any  one  of  the  four  elements  of  the  reflection  matrix.  The  elements  of 
the  coefficients  are  then  given  by 

a()  =  Id;  <  R^  [  +  RLR)_dT(RUR  +  RLL)'/I)2 
a,  =  l  d  t  ( R  y  l  -  RLR)-d?(RUR~  RLL),,/DI 

a: =  i-(Rul  +  rlr>  +  (Rur  +  rll^  ^2 

a  >  =  I  -d^*  Rl’l  “  RLR*  +  d  1*  Rl'R  ~  R1.L*I  d)l 

where 

d  1  ~  <0 pjj  -  0c)/t  =  -  (0 -  0c),/t 
d  =  (0^  ^  -  0cl  t  -  -  (0j  |  -  0l.)  t 
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HR'S* 


D2  =  2d1d2(d^-dJ). 

The  letter  subscripts  refer  to  the  four  comers  of  the  rectangle  (figure  1 ). 

A  measure  of  the  adequacy  of  the  interpolation  series  is  needed  so  as  to  determine 
which  rectangles  need  to  be  subdivided  into  smaller  rectangles.  For  this  purpose  the  two 
eigenvalues  of  the  reflection  matrix  are  used  since  it  is  inappropriate  to  require  a  close 
relative  tolerance  on  weak  elements  of  the  matrix.  They  are  given  by 

R±=  {Rn  +R22+  [(RH  -R22)2  +  4Ri2R2lI,/2}/2  ■ 

For  purposes  of  determining  the  adequacy  of  the  interpolation  series,  the  values  of 
the  R±  found  from  the  elements  of  the  reflection  matrix  obtained  from  the  full-wave 
solution  are  compared  to  the  values  of  R±  found  from  elements  of  a  best-fit  second-order 
interpolation  series.  For  this  series  a2  is  set  to  0,  and  aj  is  found  from 

(al)LS=  [d^(RUL- RLR) +  d2(RUR  -  RLL)]/[2(d*dj  +  d*d2)]  . 

The  measure  of  adequacy  is  then  taken  to  be  inverse  to  the  largest  of  the  eight  values 

e+  =  H(r+)lS  ~  (R+>FW]  /[(r-PlS  +  (R+>FW1 1 
e_  =  l[(R_)LS  ’  <r-)fw1  /Kr+)lS  +  (R+>F\vl  I 

at  the  four  corners  of  the  rectangle. 

The  same  formulation  for  interpolation  is  used  for  the  complex  logs  of  the  magneto¬ 
ionic  reflection  coefficients.  In  the  above  equation  R  is  then  replaced  by  £n(Roe).  However, 
for  the  case  of  complex  logs  of  reflection  elements  there  is  an  ambiguity  in  cycles  of  phase 
to  be  resolved.  For  this  purpose  it  is  noted  that  this  phase  ambiguity  implies  an  ambiguity 
in  the  value  of  the  coefficient  a3  such  that  its  value  might  be  different  by 

=  (d  j  n  j  -  d n -> )27ri/D  j  , 

where  nj  and  n-i  are  integers.  Those  values  of  nj  and  n-i  are  chosen  for  which  la 3  +  Aa^l 
is  the  smallest.  Then  Cn(  R^eR)  and  Cn(R^)  are  modified  by 

Cn(RyeR)  £n(RyeR)  +  n  j  27ri 

?n(Rye^)  R +  n')^7r*  • 

There  is  also  an  ambiguity  in  a2  such  that  the  value  of  a-»  might  be  different  by 
Aa2  =  4nffi/D->  , 

where  n  is  an  integer.  That  value  of  n  is  chosen  for  which  |a->  +  Aa->|  is  smallest.  Then 
fin(RyCR)  and  8n(R^)  are  modified  by 


fin(R^fR)  «-  fin(RjjeR)  +  2nm 
fin(R  Jl)  +-  finCR^)  +  2nm 


VI.  EIGEN  ANGLE  SEARCH  USING  TRIANGULAR  MESH  UNITS 

The  search  scheme  formerly  used  in  MODESRCH  has  been  replaced  with  an 
algorithm  using  a  mesh  of  equilateral  triangles.  A  mesh  of  squares  was  used  in  the  earlier 
algorithm.  The  basic  principles  are  the  same  in  the  two  versions.  The  mesh  of  equilateral 
triangles  is  used  because  it  leads  to  simpler  coding  and  because  the  function  need  only  be 
evaluated  at  three  corners  of  each  mesh  unit  considered  rather  than  at  four. 

As  with  the  earlier  version,  each  of  one  or  more  lines  of  constant  phase  is  followed 
from  one  point  on  the  perimeter  of  the  search  area  to  its  point  of  exit,  provided  that  any 
such  line  passes  through  the  search  area.  The  lines  are  defined  by  the  condition  that 
lm(F)  =  0. 

The  chosen  rectangular  search  area  is  shown  in  figure  3  along  with  the  triangular 
mesh  and  hypothetical  lines  of  Im(F)  =  0.  The  search  begins  at  the  upper  left  comer  of  the 
mesh  pattern  and  proceeds  counterclockwise  around  the  perimeter  of  the  pattern.  A  phase 
line  of  lm(F)  =  0  is  detected  by  a  change  in  sign  of  the  lm(F).  Note  that  on  one  side  of  the 


+  -  -  .  ♦ 


Figure  3.  Search  rectangle  and  associated  triangular  mesh  pattern  shown  with  hypothetical  lines 
of  Re(F)  =  0  and  lm(F)  =  0  and  eigenangles. 


phase  line  lm(F)  >0  and  on  the  other  side  Im(F)<0.  Each  triangle  along  the  line  is  checked 
for  a  line  of  Re(F)  =  0  also  passing  through  the  triangle.  Such  a  crossing  of  phase  lines  indi¬ 
cates  the  presence  of  a  zero  of  the  function  F,  hence  a  waveguide  eigenangle.  A  Newton 
iteration  is  then  used  to  locate  the  eigenangle  in  or  near  that  triangle. 


VII.  EXAMPLE 

The  example  given  is  that  of  a  nighttime  ionosphere  model  at  a  frequency  of  60  kHz. 
In  terms  of  parameters  used  in  MODESRCH  the  electron  density  is  specified  as  /3  =  0.5, 
h'  =  87.  The  codip  angle  is  30  deg,  and  the  geomagnetic  azimuth  of  propagation  is  35  deg 
east  of  north.  The  full-wave  integration  was  carried  out  over  a  height  range  of  97  to  77  km 
altitude. 

The  elements  of  ( R  +  1 )/C  are  listed  in  table  1 .  These  are  the  variables  that  were  used 
in  the  previous  version  of  MODESRCH  for  all  search  rectangles  in  the  interpolation  required 
for  approximate  starting  values  of  eigenangles.  They  are  used  in  the  present  version  only  for 


Table  1.  Reflection  coefficients. 


0 

88 

86 

84 

82 

80 

78 

76 

74 

72 

70 


(nR11  +  l)/C 

11.887+49.143 
1 1.280H5.771 
9.948+43.115 
8.488+il  .367 
7.197+iO.l  1 1 
5.684-i0.900 
3.958-il  .058 
2.859-i0.383 
2.703+i0.305 
2.921+i0.355 


11R1 


,/c 


-2.274-H9.323 

1.046+45.148 

2.505+12.021 


1 .981  -iO.525 


036-il  .258 
-  805+i0.065 
.1 72+i0.684 
.717-10.141 
-,038-i0.768 
-.771+10.082 


68 


2.790+i0.065 


-  200+i0.697 


lRll/< 
-.634+17.358 
1.498+i3.861 
2.353+41.237 
1 .631  -iO. 847 
-.277  -i  1 .266 
-.910+40.1  50 
,168+i0.759 


(1Ri  +  1  )/C 

13.380-il  5.865 
6.028 -i7.086 
2.216-40.907 
3.3  1  8+43.994 
7.848+13.892 
8.267-11 .395 
3.3  8  5-43.183 


.764-10.125 
-0.3 2- i0.8 1  7 
-  839-i0.072 
-.1  79+i0.802 


.754+10.037 
2.858+i2.483 
5.094+10.467 
3.182-11  95: 


66 

64 

62 

60 

58 

56 

54 


2. 4  5  2+ iO. 007 
2 .3 1 7+i0.096 
2.302-i0.014 
2  049-i0.222 
1 .665-10.139 
1 .569+40.167 
1.782+i0.263 
1  890+i0.009 
1 .64  7- 40  2  3  8 


,573+i0.292 
,336-i0.447 
-,345-i0.336 
-.3 1 5+40  265 
,200+i0.292 
,275-i0.150 
-.1 1  l-i0.265 
-.262+10.072 


,729+i0.284 
.373-40.652 
-,579-i0.438 
-  490+j0  503 


.415+10.545 

.608-10.312 

-,188-i0.668 


.714-10.653 
1.491+41.521 
3.454+i0.856 
2.903-1 1  150 
.940-10.978 
.74  7+i0.707 
2.169+il  .053 
2.649-i0.24‘7 


52 

50 


,023+i0.259 


705+10.032 
1  56+10.699 


1.530-i0.969 
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some  selected  rectangles  and  for  very  low  ionosphere  models.  In  table  2  are  listed  the  values 
of  the  natural  logs  of  the  elements  of  the  magnetoionic  reflection  coefficients.  These  vari¬ 
ables  are  used  for  the  interpolation  used  in  the  newer  version  of  MODESRCH  for  most 
search  rectangles. 

Table  3  shows  the  magnitudes  and  phase  values  of  the  magnetoionic  reflection 
matrices  at  “given”  points  on  the  real  0-axis.  These  points  are  at  2-deg  intervals  from  50  to 
88  deg.  Note  that  leRel  (upper  left  element)  is  much  larger  than  leRQl  (lower  left  element) 
for  all  values  of  0  and  that  l0R0l  (lower  right  element)  is  larger  than  l0Re I  for  the  stronger 
modes.  This  indicates  that  reflection  is  principally  from  two  distinct  complex  heights. 

These  effective  heights  are  also  listed  in  table  3. 

Propagation  parameters  from  the  previous  version  of  MODESRCH  for  the  highest 
44  modes  are  listed  in  table  4.  In  table  5  are  given  the  propagation  parameters  from  the 
newer  version  of  MODESRCH.  The  eigenangles  agree  to  ±0.001  deg.  The  agreement  for 
attenuation,  phase  velocity,  excitation  factor,  and  polarization  is  also  more  than  adequate. 


Table  2.  Complex  logs  of  magnetoionic  reflection  coefficients. 


0 

eRe 

eRo 

oRe 

oRo 

88 

-.215+00.182 

-2.743+127.490 

-2.01  7+i  1 7.1  19 

-.554+120.463 

86 

-.149+120.487 

-2.71  l+i26.831 

-2.063+116.672 

-.656+120.178  J 

84 

-.1 73+128.525 

-2.651+i25.860 

-2.204+116.050 

-.806+119.765 

82 

-,147i+27.399 

-2.546+i24.69 4 

-2.498+il5.333 

-.992+H9.259 

80 

-.128+126  156 

-2.418+i23.421 

-3.046+114.624 

-1.206+118.677 

78 

-.1  14+124.822 

-2.303+i22  068 

-3.966+114.426 

-1.436+118.021 

76 

-.097+123.425 

-2.195+120  63" 

-3  549+1 1 4.726 

-1.659+117.288  i 

74 

-.075+121  995 

-2  061+119  |6  | 

-3  103+113.892 

-  1  845+1 1 6.477 

-  *> 

-.050+120.549 

-1  908+1  p  (.89 

-2  84’+i  12  816 

-1 .475+,  15 .61 1 

70 

-.027+119.091 

- 1  ,-66+i  1  <>.226 

-2  '09+,l  1  '04 

-2.063+114.724  ; 

68 

-  002+1  r  626 

-1  639+114  ''4 

-2  683+1 1 0  642 

-2  152+113.841  I 

66 

031+116  161 

-1  5()2+i  t  3.2'U 

-2  8 1 8+|9  (,58 

-2.286+M  2.457  | 

64 

.073+;  14.7  13 

- 1  346+i  1  1  804 

-3  148+18.8  18 

-2.490+112.067  ; 

i 

62 

.123+113  295 

-1. 184+110. 376 

-3  524+iX  327 

-2.747+1 11.162 

60 

.176+11 1.421 

-1  .032+18  989 

-3.342+P  971 

-3.01 5+110.235  j 

s.x 

.21  6+1 1  0  60' 

-O  oi)S  +  r  i,4  s 

-2  "68+ 1'  214 

( 

-3  1 9S+|9  421 

56 

.203+10.350 

-  8 24 ♦  j6  332 

-2. ('24  +  i(i  3  1  ' 

-3  236+18.44'  i 

54 

.130+18.1  10 

-783+p  or 

-2.456+i5.423 

-3.187+17.584  ; 

52 

.027+16.868 

-.7*2+13.69? 

-2.468+14.519 

-3  1  24+16. 67'’  i 

50 

-  O80+1 5  (02 

-  ~  1 8+i."  is ; 

-  2  ^  '1,+  1  '  M\S 

-3  100+0  69(1  | 

' 

Table  3.  Magnitude  and  phase  of  magnetoionic  reflection  coefficients  and 
apparent  complex  heights  of  reflection. 


/.863  .082 

1.078  .371 


/  129.9°  158.5° 

l  -25.2°  23.5° 


/.880 
l  .089 


/  58.6°  117.9° 

(  -98.1°  -9.9° 


/.893  .019 

(.100  .238 


/  -17.8°  106.5°\ 

(-175.6°  -47.5  y 


/  -97.8°  123.7° 

(  102.4°  -89.5° 


/-1 79.8°  76.0° 

(  17.8°  -135.9° 


(97.3°  14.3° 

-66.5°  174.4° 


.973 


/  13.8°  -49.4° 

(  - 1 50.3  123.6° 


998 

194 


/  -70.1°  -110.3° 

(  125.3°  73.0° 


e  =  66° 


/ 1  03  I  060  \ 

\  223  102  ) 


-154.0°  - 1 66.6°  \ 

40.3°  22.4°  ) 


82.80+i0.1 1 
79.32 -i0.83 

88.74+i0.25 
8 1.07-11.47 

89.08+i0.31 
82.31 -i  1.94 

90.68+i0.26 

83.27-i2.32 

91 ,96+iO.l  8 
84.l6-i2.60 

92.96+10.17 

85.08-12.69 

93. 64+  iO.  23 
86 . 09 -i  2. 46 

94.06+i0.28 
87.00-1 1.88 

94.39+i0.29 

87.56-il.24 

94.73+10.28 

87.76-i0.98 

95.03+i0.35 

87.85-11.30 

95.19+10.47 

88.03-12.10 


Table  3.  Continued. 
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/  123.0° 
y  -43.7° 


41.7° 

-125.5° 


-112.3 

78.0° 


96.7 
—1 33.6° 


-101.1 
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Table  4.  Propagation  parameters  using  full-wave  solution  values  of  R  for  each  eigenangle. 
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VIII.  SUMMARY 


A  formulation  has  been  described  that  has  made  it  possible  to  find  LF  earth-iono¬ 
sphere  waveguide  mode  constants  at  a  much  reduced  cost  and  very  low  frequency  (VLF) 
mode  constants  at  less  cost  than  previously.  The  most  relevant  cases  are  those  for  which  a 
large  number  of  modes  are  required,  such  as  for  short  distances  from  the  transmitter,  night¬ 
time  propagation,  or  elevated  antenna. 

In  contrast  to  the  previous  version  of  MODF.SRCH,  the  mode  constants  are  func¬ 
tions  of  interpolated  values  of  ionosphere  reflection  coefficients  rather  than  of  reflection 
coefficients  found  by  a  full-wave  solution  for  the  eigenangle  associated  with  each  mode. 

The  given  values  for  the  interpolation  are  found  with  the  full-wave  solution  only  at  selected 
points  in  the  0-plane.  Examples  indicate  that  the  interpolated  values  are  more  than  ade¬ 
quate  for  finding  the  values  of  the  propagation  mode  constants. 

A  formulation  is  presented  that  appears  to  be  effective  in  separating  waves  reflected 
from  two  rather  distinct  complex  heights  in  the  ionosphere.  The  advantage  is  that  compli¬ 
cations  from  the  interference  of  the  two  waves  tend  to  be  avoided.  Values  of  magnetoionic 
reflection  coefficients  found  in  the  example  indicate  that  adequate  separation  of  compon¬ 
ents  was  accomplished  for  all  values  of  the  incidence  angle.  0.  required.  The  separation  need 
not  be  complete  for  successful  interpolation.  It  was  more  than  adequate  in  the  example, 
however,  and  especially  complete  for  values  o!'0  near  90  deg.  where  eigenanglex  sometimes 
are  difficult  to  find  using  the  previous  version  of  MODESRCH. 

The  separation  cannot  be  carried  out  in  the  vicinity  of  branch  points,  where  mag¬ 
netoionic  engenvectors  cannot  be  clearly  distinguished.  In  the  example,  no  branch  points 
occurred  in  the  vicinity  of  the  eigenangles.  The  formulation  does  allow  for  these  branch 
points  to  be  located,  however,  and  their  occurrence  should  present  no  problem  since  the 
procedure  allows  for  defaulting  to  a  method  of  solution  more  nearly  like  that  used  in  the 
previous  version  of  MODF.SRCH  in  the  neighborhood  of  a  branch  point.  Such  cases  are 
more  likely  to  occur  for  cases  of  daytime  ionosphere  models  than  for  nighttime. 

The  example  is  for  a  nighttime  ionosphere  at  oO  kll/  and  for  propagation  over  sea¬ 
water.  Forty-four  eigenangles  and  associated  mode  constants  were  found  in  the  interval  of 
50  to  90  deg  in  the  real  part  of  the  incidence  angle.  0.  1  he  cost  was  onlv  about  one-sixth 
that  required  for  finding  the  same  eigenangles  and  mode  constants  using  the  previous 
version  of  MODESRCH  Eigenangles  differed  at  most  In  0.001  deg. 
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